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Abstract 

A scheme for solving quasiclassical string equations is devefopped to prove that genus-zero Whitham 
hierarchies describe the deformations of planar domains determined by rational conformal maps. This 
property is applied in normal matrix models to show that deformations of simply-connected supports of 
eigenvalues under changes of coupling constants are governed by genus-zero Whitham hierarchies. 
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1 Introduction 



Conformal mapping methods have been effectively apphed in the analysis of interfacial free-boundary prob- 
lems involving planar domains [T. They have provided many exact solutions [H-|H which stimulated the 
research on possible underlying integrable structures. Thus, Wiegmann and Zabrodin discovered 0-10] that 
deformations of simply-connected domains with respect to changes of their exterior harmonic moments, 
treated as independent variables, are described by the dispersionless Toda hierarchy. They also formulated 
an algebro-geometric analysis (7j-|Hl of the deformations of multiply-connected domains in terms of Whitham 
equations for Abelian differentials. 

Recent research has shown [5]-|lUj that many exact solutions of Laplacian growth models correspond 
to a special type of domains called algebraic or quadrature domains. In the simply-connected case the 
complement of a quadrature domain D is the image of the exterior of the unit disk under a conformal map 
given by a rational function 

A^o k Ns 

n=0 s=l n=l ^ 

where the coefficient r is a positive number and the k poles as ^ lie inside the unit circle. In this work we 
prove that the deformations of rational conformal maps under changes of the parameters {r,uo^n,Us,m,CLs) 
such that z{aj^) are kept constant, turn out to be described by a solution of the genus-zero Whitham 
hierarchy W(n) with n = 2k + 2 punctures It should be noted that according to a recent general result 
by Takasaky ^2 the Whitham hierarchy W(n) is the quasiclassical limit of the n-component KP hierarchy. 

Our analysis is based on solving a system of quasiclassical string equations which leads to the character- 
ization of the conformal map as a function z(w, t) of the Whitham times t. For k = our result agrees 
with [3] (see also |13J) since W(2) is the dispersionless Toda hierarchy. However, for k > 1 the analysis of 
[H] does not apply because the Schwarz function of the boundary of D has poles outside D, consequently 
there are infinite exterior harmonic moments of D different from zero and, since z{w) depends on a finite 
number of parameters only, these harmonic moments are not independent variables. 

Deformations of quadrature domains naturally arise in the analysis of partition functions oi N x N 
normal matrix models 

Zj^= [ e^^(*^>^')dMdMt, (2) 



with quasiharmonic potentials W{z, z) := —zz + V{z) + V{z) of the form 

No+l k Ns-1 

v{z) = y: ^"io,n + E ( - *^.o iog(. -ps)+y: 77^) • 

n=l s=l n=l ^ 

Integrating over eigenvalues and ignoring normalization factors, the partition function reduces to 

jj\\z^-z,\^e^^^'^^^^-^^^\{d^z,, (4) 

i>j 3 



In the large N limit {N oo, hN fixed) the eigenvalues densely occupy a bounded quadrature domain D in 
the complex plane {the support of eigenvalues). As a consequence of our analysis we prove that for simply- 
connected supports of eigenvalues the corresponding rational conformal map z = z{w, t) as a function of the 
coupling constants t of the partition function represents a solution of the Whitham hierarchy W(2A; -|- 2). 



2 



Figure 1: Right-exteriors of 7^ and 



2 String equations in Whitham hierarchies 

The elements of the phase space for a genus-zero Whitham hierarchy W(M + 1) are characterized by M + 1 
punctures Qa, (a = 0, . . . ,M), where go := oO; of the extended complex p-plane and an associated set of 
local coordinates of the form 



00 , 00 



zo=P 



+E^' ^i = ^.+T.dUp-<i^r. i = i,...,M. (5) 



n=l P P - n=0 



In what follows Greek and Latin suffixes will label indices of the sets {0, . . . , M} and {1, . . . , M}, respectively. 
We will henceforth suppose that there exist positively oriented closed curves in the complex planes of 
the variables such that each function z^{p) determines a conformal map of the right-exterior of a circle 
7^ := z~^{Tfj_) on the exterior of F^ (we will assume that the circle 70 encircles all the 7^) (see figure 1). 

The flows of the Whitham hierarchy can be formulated as the following infinite system of quasiclassical 
Lax equations 

dz 

= {nf,n,Za}, (6) 

associated to the series of time parameters {to,n ■ n > 1 ]ti^ri '■ i = 1,...,M, n > 0}. Here the Poisson 
bracket is defined as {F, G} := dpF dxG — dxF dpG, x := toi and the Hamiltonian functions are 

^f,n ■■= {z^){i,,+), (n > 1); Qio := - log{p-qi), (7) 



00 





where (•)(i,+) and (■)(o,+) stand for the projectors on the subspaces generated by {{p—qi)~^}^=i and {p"},^ 
in the corresponding spaces of Laurent series. These hierarchies |12| are the dispersionless limits of the multi- 
component KP hierarchies. In particular, for M = and M = 1 they represent the dispersionless versions 
of the KP and Toda hierarchies, respectively. 

In our analysis we will use an extended Lax formalism with Orlov functions 



• „^ .... . , 

n=l n>2 ^' 
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which verify the same Lcix equations as the variables Zq,^ and such, that {-Sq:, nia} = 1. The parameter 
too in (© is defined by 

M 

too '■= — ^^iO, (9) 

i=l 

The Whitham hierarchy can be formulated as the fohowing system of equations 

dza A drua = duj, Va, (10) 

where uj is the one-form 



:=^n^„dt^n- (11) 



UJ 

To see how to get from the system (jlOj) to the Whitham hierarchy, note that by identifying the coefficients 
of dp A dt^n and dx A dt^n in (ITU)) we obtain 



dza drua dnia dza dVt^n dza drua dnia dza dVt^n 



(12) 



dp dt^n dp dt^n dp dx dt^n dx dt^n dx 

and, in particular, since iloi = for (/i,n) = (0, 1), the system ((T^ implies {za,ma} = 1. Thus, using this 
fact and solving (|T^ for dt^„Za and dt^„'ma, the Lax equations for {za,ma) follow. 

A natural form of characterizing solutions of Whitham hierarchies is provided by systems of string 
equations 

(PM,m) = Po{zo,mo), ,^i^2,...,M, (13) 
[Qi{zi,mi) = Qo{zo,mo), 

where {Pa,Qa}aL(} satisfy {Pa{p,x),Qa{p,x)} = 1. Given a solution {za{p,t),'ma{p,t)) of a system (|TT?|) . if 
we denote 

Va{p,t) := Pa{Zaip,t),ma{p,t)), Qa{p,t) ■= QaiZa{P,t) , rUaip, t)) , 

it is clear that 

dVa A dQa = dVp A dQp, Va, /3. (14) 

On the other hand {Pa{p-,x)^Qa{p-,x)} = {zq,?tIq,} = 1, so that solutions of a system of string equations 
verify 

dVa A dQa = dzp A dmp, Va, (3. (15) 



Theorem Let {za{p,t),ma{p,t)) be a solution of (|13j) which admits expansions of the form (jSI) , (jHI) > © o-i^-d 
such that the coefficients of the the two-forms ()15|) are meromorphic functions of the complex variable p with 
finite poles at {qi, . . . ,qni} only. Then {zo,{p,t),mo,{p,t)) is a solution of the Whitham hierarchy. 

Proof. In view of the hypothesis of the theorem the coefficients of the two-forms ()15() with respect to the 
basis 

{dpAdtan, dtan^dti3m} 

are determined by their principal parts at g^, (^ = 0, . . . , M), so that by taking (|15j) into account we may 
write 

M 

dZa A drUa = '^{dz^ A dm^)(^^+), Va. 
fi=0 
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Moreover the terms in these decompositions can be found by using the expansions (jH)) of the functions 
as follows 

dz, A dm, = d.^ A ( f; nz]i-'dt,n + ^ + E 

n=l n>2 
1 dVn 



= d( E + log z,dt,o - E ^t^ttS) > 

n=l n>2 ^ 

SO that 

CO 

(dz^ A dm^)(^^4.) = d(^ E(^M)(M,+)dV - (1 - ^f^o) log(p - g/i)dt^o) = d( E ^/^dt^™) • 



n=l 

Thus we find 



dza A dm^ = do; = df E ^fmdtf^n), Va 



Ml™ 

and, consequently, this proves that the functions {za{p,t),ma{p,t)) determine a solution of the Whitham 
hierarchy. □ 

3 Integrable dynamics of quadrature domains 

Let us consider a rational conformal map z^w) of the form with k poles inside the unit circle and 
define 

A^O k Ns - n 

n=0 s=l n=l 

In order to establish the connection between the deformations of the conformal map z{w) and the genus-zero 
Whitham hierarchies we introduce the change of variable 

p = Rw := r w + uqq, 

where r and uoo are the first coefficients of z{w) in As a function of the new variable p the conformal 
map is normalized at infinity 

z{p) = p + 0{l/p), p^oo. 

Moreover, if we define 

qo := oo, Qs := Rbs, Qs+k ■= Rag, q2k+i ■= uoo, 

it is clear that z{p) and z{p) become rational functions of p with poles at {qo,qs,q2k+i) and {qo, qg+k-, Q2k+i) 
respectively. 

We are going to prove that deformations of z{w) with respect to the coefficients 

U := {r,Uo,n,Us,m,as), 
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such that Ps = ^{'^)\bs ^1*6 kept constant, are described by the Whitham hierarchy W(2A; + 2). To this end 
we introduce Whitham variables (^^q,, rua) on the 2k + 2 punctures of the p-plane 



zq = z, niQ = z, (near qq = 00) 



■nis = -{z - (3sf z, (near g's) 



(17) 



Zs+k 



Z-Ps' 



nis+k = {z- Psf z (near Qs+k)- 



^Z2k+i = z, m2k+i = -z, [neai q2k+i) 

It is clear that these variables are rational functions of p with possible poles at the punctures qa only. 
Moreover, they satisfy the system of string equations 

r 1 



\- Ps = rUs+kZs+k = -'m2k+i = zo 

■2c 



(18) 



-rusz'^ 



Zs+k 



+ Ps = Z2k+i = mo. 



Obviously the functions Za are of the form ©. On the other hand due to and ((TB|) it follows that the 
functions nia defined in (|T7]) verify expansions of the form (jH)) 



E n-l , too , 
ntonZQ H \ 

Zo 



n=l 



-Zs "^Z 



n=l 



^ Zs 



W — > CO, 



(19) 



rris+k - Zg^i^z 



m2k+i = -z 



E n-l , ts+kO , 
ntsj^knZ^_^_f. + \ w ^ as, 



n=l 

No+1 



E, n-l , hk+10 , 
nt2k+lnZ2k+l + + 



n=l 



Z2k+1 



where the time parameters t := (ta^n) are rational functions in u 

tan = Qan{u). (20) 

The functions Qom{u) satisfy certain constraints which can be characterized by considering the map 

C:f^Cf, Cfiw):=J{^. (21) 
Observe that in terms of the variable p 



Cf{p) = f{Ip), 



(22) 
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where Tp = r"^ /{p — uqo) + ^oo is the inversion with respect to the circle \p — uqoP = r'^- From (|17|) it is 
clear that 

Z2k+i = Czo, m2k+i = -Cmo, Zg+k = Czs, tUs+k = -Cms, (23) 

which implies 

Q2k+l,n{u) = -Qo,n{u), Qs+k,n = -Qs,n{u). (24) 

Furthermore, we can prove that 

J2Qao{u) = 0. (25) 

a 

Indeed from (jlZf) we deduce 

mo dzQ = rus dzs = z dz, m2k+idz2k+i = rris+kdzs+k = -z dz. 

Hence 



27ri ^ Qao = X] f '^c.dza = f z dpzdp = 0, 



where we have taken into account that z dpZ is a rational function of p with poles at the punctures qa only, 
and the fact that 

^7„~0 in C\{gi,...,g2fc+i}- 

a 

Notice that due to ()24() the constraint ()25() can be rewritten as 

ImQoo + ^ImQsoiu) (26) 

s 

Under appropriate conditions one can determine u as a function of {t,Ps,Ps)- To this end we consider the 
system 

tan — Qan{'^^: 

(27) 

^Ps = z{bs,u), 
where the time parameters 

t:={ta,n), a = 0,...,2k + l;n = 0,...Na; Nq = N2k+i = No + I, = N^+k = Ns - 1, (28) 
are assumed to satisfy 

t2k+l,n = —to,n, ts-\-k,n = —ts,n, ^aO = 0. (29) 

a 

Firstly, we observe that u constitutes a set of 

2{No + Y,Ns + k)+3, (30) 

s 

real variables given by r and the real and imaginary parts of {uQ^n,Us,rmCLi)- On the other hand, in view 
of H24|) and H25() we may ignore the equations corresponding to a = s + /c, 2fe + 1 in ()27() . In this way the 
system (|27|) reduces to 2(A'^o + 2 + J2s + k) real equations, but due to (|26jl one of them is a consequence 
of the others. Therefore, we are lead to a system of equal number of equations and unknowns which under 
appropriate conditions will determine u, and consequently (za,iTLa) as functions of (t,(3s). 
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In this way we have determined a rational solution {za{p) , maip)) of the system of string equations 1)18(1 
which depends on (t, . . . , and satisfies the asymptotic conditions © and Q. Therefore, from the 
above Theorem we conclude that this solution evolves with respect to t according to the Whitham hierar 
chy W(2/c + 2). 

It is interesting to notice the following identity involving the times tao and the area T of the domain D. 
Let r be the positively oriented boundary of D and take small positively oriented closed curves Fg around 
the points of the z-plane, then we have that 




(31) 



4 Exact solutions 

The above analysis can be used to generate solutions of W(2k + 2) of the form (0). Let us consider the case 
in which only simple poles arise 

k k _ 

X - Vs _ r _ \ - VgW 

z = rw ' ^- ' - . i 



, z = — + UQ + y — 

w — as w 1 — ttsW 

By identifying the coefficients of and Z(^^ in the expansion of mo = z = z ^ oa (w ^ oo), we get 

k _ I k _ \ 

and identifying the coefficient of z'^'^ in the expansion of nis = —{z — Ps)'^z as — > oo {w bs) yields 

Finally, the equations z{w)\uj=b^ = (3s read 

a, 1 — a, a,/ 



-,2 



k _ 
Vs'CLs 



The system (|32() - (|34() determines {r,uo,Vs, as) in terms of (to,0) *o,i; ^5,0; /^s)- 
Examples 

A solution of W(4) is obtained by deforming the conformal map (aircraft wind) 



V r _ vw 
z = rw + uo-\ , z = \-uo + 



w — a w 1 — aw 
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Figure 2: Solution corresponding to k = 1 



The solution of the corresponding equations (jH^ - lj^H) is given by: 



- AtoA 



r{P-to,i){l-A) 



4 2 



Uo 



toflA + rHo^i {(3 - to,i)(l -A) + M^A - to,( 



with 



r2 

2|to,oP^' 



1/3 - to,i P (2 A3 - 3 + ^) + (to,o + ti,o)(^ - 1) + (tofi + to,o) A ' 



where A = \a\ is implicitly defined by 



2|/3-to,i|^^^+ 2(to,o + ti,o)-2(to,o+to,o)-5|/3-to,i|' |/3-to,i|'^' 



|2 i4 



+4 (to,o + to,o) - (to,o + ii,o) + 1/5 - to.iP 1/3 - io,iP^^ 



+ (to,o + h,of - (to,o + to,o? + 4|to,oP - 2(to,o + to,o)\P - to,: 



+2(to,o + hfiW - to,i? to,if\A^ - 2(^0,0 + h,ofA + (to,o + ti,of = 0. 

It can be proved that a reduction of the system H32() - H34() is obtained by setting 



z = rw 



Vs ^ Vs 

w — ac, w — 



w 



VgW ^ VsW 



s=l 



a^w 1 — 



and by assuming 



too, toi G IK, tg+i o — ts,0, Ps+l — (3s. 
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Figure 3: Solution corresponding to k = 2 



In the simplest case / = 1 the reduced system of implicit equations reads 

Vl Vl ( Vl Vi\ 

u — =to,i, rlr =^0.0, 

ai ai \ oi^ af J 



r Vl Vl 

''''.^"(l-af)2"(l-aiai)2 



r , , viai viai 

— + U+ — + Z2 = A 

oi i — aiQi 1 — af 



tin, (36) 



and determines a solution of W(6). 

Figures 2 and 3 exhibit deformations of domains with respect to changes of the area T such that the 
Whitham times tso are kept constant. Observe that the evolution of the boundary develops cusp-like 
singularities ^SJ. 



5 Whitham times and coupling constants of normal matrix models 

Let us assume that the support of eigenvalues of the normal matrix model © is a simply-connected domain 
D. For example this is the case if k = 1 [S]. Then by applying the saddle point method to the large N limit 
of (lU it follows that 

T.// N 1 ff dz'Adz' 

Let us consider now the rational conformal map z{w) associated to D, from 1)3 7j) it follows that the function 
z{w) can be extended as a meromorphic function of z outside D by 

N T^// N 1 f f dz' Adz' 

In fact z{z) represents the Schwarz function of the boundary F of D. By using © and H17() one can rewrite 
(jHH|) as 

^0+1 k Ns-i j^/ . J./ 

Z{Z) = ^ nz^-\n - E( E ^^s^'tsn + Z^tso) + IT^ ,_ ■ (39) 

1 1 1 ^'^'^ J J D Z Z 

n=l s=l n=l 
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By comparing this identity with our definition H19j) of the Whitham times and by taking into account (|31() , 
we conclude that the Whitham times coincide with the coupling constants of the normal matrix model. 

Finally, we notice that by using the strategy deployed in the proof of the Theorem of Section 2 and 
by taking (|39() into account, it follows that the deformations with respect to the parameters and (is are 
described by the flows 

^ = + E'^^^.-T^ + (^^ - ^%,N.-i{{z^%,^),z^], (40) 

and 

r, O Ns-2 „ 

= ts+k,Ojr. + V nts+k,nTr. + [Ns - 'i^)ts+k,Ns-l{{Zs+k){s+k,+), Za}- (41) 

OPs Ols+k,l Ols+k,n+l 
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